Velocity-force characteristics of a driven interface in a disordered medium 
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Using a dynamic functional renormalization group treatment of driven elastic interfaces in a 
disordered medium, we investigate several aspects of the creep-type motion induced by external 
forces below the depinning threshold f c : i) We show that in the experimentally important regime of 
forces slightly below f c the velocity obeys an Arrhenius-type law v ~ exp[— U(f)/T] with an effective 
energy barrier U(f) oc (f c — /) vanishing linearly when / approaches the threshold / c . ii) Thermal 
fluctuations soften the pinning landscape at high temperatures. Determining the corresponding 
velocity-force characteristics at low driving forces for internal dimensions d = 1,2 (strings and 
interfaces) we find a particular non-Arrhenius type creep v ~ exp[—(f c (T)/f) >1 ] involving the reduced 
threshold force f c (T) alone. For d = 3 we obtain a similar v—f characteristic which is, however, 
non-universal and depends explicitly on the microscopic cutoff. 

PACS numbers: 05.20.-y, 64.60.Cn, 74.60.Ge, 82.65.Dp 



I. INTRODUCTION 



The influence of disorder on the static and dynamic 
properties oLelastic systems has been intensely studied in 
recent yearscm. Various physical systems including 
lattices in superconductors^, domain walls in magnetsc 
and charge density waves in solidaj significantly change 
their properties upon introducing even a small amount of 
disorder. Subject to a disorder landscape, these systems 
transform to a glassy state characterized by a nontrivial 
scaling of the displacement correlation functional ajad-a 
vanishing linear response to external driving forcesaliil, 
e.g., the current induced Lorentz force acting on vortices 
or the magnetic field driving the domain walls in mag- 
nets. The determination of the velocity-force character- 
istics of a driven elastic manifold subject to a disorder 
landscape is a challenging problem: while the behavior 
at small distances and large drives is amenable to pertur- 
bation theory, the most interesting long distance/weak 
drive regime can only be attacked via non-perturbative 
methods. In this paper, we consider some aspects of the 
creep-type dissipative motion of a driven elastic interface 
with d internal dimensions, moving along one transverse 
direction in a disorder landscape (d + 1-dimensional ran- 
dom manifold problem). 

Depending on the value of the temperature T and the 
external force / several regimes can be distinguished (see 
Fig. |l|): At T = 0, the velocity v is zero as long as / does 
not exceed the critical force f c , whereas for / > f c the 
system starts moving, v(f) ^ 0. In particular, one finds 
v(f ) oc (/ — fc) 13 near the threshold (the depimung tran- 
sition), with a nontrivial critical exponent For 
large drives / 3> f c the disorder becomes irrelevant and 
the velocity- force characteristic turns linear, v ~ f/rj, 
with rj the friction coefficient characteristic of the dissi- 
pative dynamics. 



At finite temperatures T > 0, thermal fluctuations in- 
duce a creep-type motion resulting in an exponentially 
small but finite velocity even below threshold f < f c (see 
Fig. 0). At small drives / — > an Arrhenius-type law 
v(f) oc exp {— U(f)/T} holds, with a diverging activa- 
tion barrier [/(/—> 0) — ► oo (glassy response). 
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FIG. 1. Velocity- force characteristic of a driven interface. 
The thick solid line is the zero-temperature result with a 
threshold force f c below which the velocity vanishes. Beyond 
f c the velocity first rises following the scaling law v oc (/ — fc) 13 
and then crosses over to the linear dissipative regime with 
v oc /. The thin line shows the behavior at finite but low 
temperatures with a creep regime at low forces / <C fc, 
v oc exp[— (f/ c /T)(/ c //) M ]. Close to threshold we find that the 
creep barriers scale linearly in /, v oc exp[— (U c /T)(l — f/fc)]- 
At high temperatures (dash-dotted line) thermal fluctuations 
become particularly important in dimensions d — 1,2; the 
threshold force f c (T) is strongly reduced by thermal fluctu- 
ations and we find a non-Arrhenius glassy response at small 
drive with v oc exp[— (/ C (T)//) M ] determined by the renor- 
malized critical force alone. 
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Close to threshold / ~ f c one may distinguish two in- 
teresting regimes: i) fixing the force / at its critical value 
f = f c , thermal fluctuations smooth the depinning tran- 
sition and the velocity vfML^-pt T 1 / 5 is expected to scale 
as a power of temperaturer 6 H 17 l. ii) fixing the temperature 
T and increasing the force towards threshold / — > f c , a 
creep-type response is expected with a vanishing activa- 
tion barrier U(f — > f c ) oc (f c — f) a — » 0; here we are 
interested in the second situation-. 

The scaling theory of creepQlHl predicts that U(f — > 
0) cx C/ C (/ C //) M , with a characteristic energy scale U c set 
by the disorder landscape. On the other hand, when the 
force / approaches f c from below one expects that the 
barrier behaves like 



t/(/)oe(/ c -/r, 



(1) 



with a an exponent depending on the dimensionality of 
the space and the elastic manifold. The parameter a de- 
termines the relaxation of magnetization in superconduc- 
tors at currents close to the critical one. In his original 
description of magnetic relaxation, Andersonta assumed 
that a = 1 and explained the famous logarithmic decay 
of the magnetic field trapped inside a superconductorliia. 
Note that the regime / ~ f c is an experimentally im- 
portant one: the observation of the system response at 
small driving forces / -C f c involves long relaxation times 
oc exp[(J7 c /T)(/ c //)' i ] and, hence, this regime is more dif- 
ficult to access experimentally. More quantitatively, the 
maximal creep barrier U that can be observed after a 
waiting time t is given by U(t) « Tln(l +t/to) with to a 
characteristic time scale involving details of the critical 
statec3u. For the vortex creep problem this time scale 
typically is of order 10~ 6 s and thus the experimentally 
attainable value of U/T is limited to factors ~ 30-40. 

From a theoretical point of view the calculation of the 
barrier exponent a in (Q) near criticality still remains a 
problem. In fact, one may expect that critical fluctua- 
tions of the manifold near the threshold will affect the 
creep motion. In this paper we study the behavior of 
U(f) near the threshold using dynamical renormaliza- 
tion group theoryta and show that if the pinning of the 
manifold is due to a short-range correlated random po- 
tential (e.g., due to point-like impurities) the effective 
barrier behaves as 



U(f) *U C (1- f/fo) 



/-/c 



(2) 



with U c a characteristic energy scale set by the disorder 
landscape. This result is independent of the dimensional- 
ity of the maaifold and confirms the original assumption 
of AndersonES. 

In addition, we investigate creep at high tempera- 
tures, again using the dynamical functional renormal- 
ization group technique. In this case, the dimensional- 
ity of the manifold is particularly important: It is well- 
known that the mean thermal displacement (« 2 }th of a 
manifold with internal dimension d > 3 is bounded, the 
maximum displacement depending on the microscopic 



short-scale cut-off of the elastic system. Strings and sur- 
faces (d = 1, 2), however, exhibit thermal fluctuations 
([u(z)-u(0)] 2 } which grow unboundedly with separation 
z. At high temperatures, the manifold probes an effective 
disorder landscape averaged over thermal displacements 
which are only bounded through the disorder-induced 
pinning at large scales, resulting in a strongly reduced 
disorder strength. In particular, the critical force f c (T) 
is found to decrease as a power law with increasing tem- 
perature, f c {T) ~ f c {T dp /T) K with k = 7 (2) in d = 1 (2) 
dimensions (see Fig. [I]). The characteristic temperature 
Tdp determining the crossover from the low to the high 
temperature regime is given by the bare disorder energy 
scale U c , Tdp = U c . For the effective barrier in d = 1, 2 
we find 



tf(/)~T[(/ c (T)//)' 4 -l] 



(3) 



depending only on the renormalized pinning force f c {T), 
confirming the results obtained previously via scaling es- 
timates, see Ref. ||. Similarly to the low T creep, the 
exponent fj, = (d + 2£ — 2)/(2 — £) again is determined 
through the static roughness exponent £. The large ther- 
mal fluctuations modify the characteristic energy scale 
U c of the problem to U c — > T, leading to a peculiar non- 
Arrhenius form for the v — /-characteristic in the high T 
creep regime. For d = 3 the velocity-force characteristic 
also takes a non-Arrhcnius form but with an exponent ad- 
ditionally modified through the temperature dependence 
of the creep barrier U C (T) 3> T. Note, however, that 
Eqs. (||) and (||) make sense only if the temperature T is 
still small enough to produce an exponent C/(/)/T»l. 

In Sec. II below we will first analyze creep near the 
critical force and derive Eq. (0) while section [II is de- 



voted to the study of creep at high temperatures with a 
derivation of the result Eq. (0) for the creep barrier. 



II. CREEP NEAR THRESHOLD 

In the following we will concentrate on the case of 
(d+ l)-dimensional elastic media, with d internal dimen- 
sions and one single transverse direction. Typical realiza- 
tions are strings confined to a plane (a (l + l)-dimensional 
manifold) or two-dimensional membranes embedded in 
three-dimensional space (a (2 + l)-dimcnsional manifold). 
These models describe domain walls in thin film and bulk 
random magnets, for example. The motion of the elastic 
manifold is governed by the equation 



■qd t u = csj 2 z u + f pin (u, z) + C(z, t) + f, 



(4) 



where the friction and external driving forces are given 
by r/dtu and / respectively, and the additional forces act- 
ing on the manifold are those due to elasticity cV^w, 
pinning f p - m (u, z), and thermal fluctuations £(z,t); rj 
and c denote the viscosity and the elasticity per unit 
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volume. We assume that the pinning force is a Gaus- 
sian random variable with zero mean and a correlator 
(f(u,z)f(u',z')} = A{u - u')6 d {z - z') of width £, the 
typical length scale of the disorder landscape. The statis- 
tics of the stochastic force £(z, t) is Gaussian as well and 
the correlator is related to the viscosity r] and the tem- 
perature T via (C(z,i)C(z',i')) = 2r 1 TS d (z - z')S(t - t'). 

The calculation of the average velocity v = (dtu) 
as a function of / and T is a difficult problem in the 
creep regime / < f c since most of the time the mani- 
fold is pinned by the random potential and only rarely 
a strong thermal fluctuation will drive it into a neigh- 
boring metastable state. Obviously, this type of motion 
cannot be described perturbatively. However, it can rig- 
orously-be proven that the velocity-force characteristic is 
uniquecil. 

A powerful method to study random elastic-manifolds 
is the functional renormalization group (FRG)E3 _ Ei.with 
various extensions dealing with finite temperature^ and 
velocityE£rli3. For dimensionalities d of the manifold 
larger than four, the effect of disorder can be taken into 
account perturbatively, whereas in less than four dimen- 
sions, an e-expansion allows to study the properties of 
the system at small e = 4 — d. The FRG has provided 
numerous results in the investigation of static and dy- 
namic propertieg-pf elastic manifolds: The static wa&* 
dering exponentE 2 ] £ as well as the dynamic exponents 
z have been determined for different types of disorder. 
Furthermore, the depinning transition at T = has been 
analyzed and the critical exponent-^ ia-the depinning law 
v cx (/ — f c )P has been calculatedl!3ii3. 

TJic dynamical extension of the FRG by Chauve et 
al.tB allows to investigate the creep regime and confirms 
the creep law U(f) cx U c (f c / /) M derived earlier via scal- 
ing arguments. In addition, it turns out that the method 
allows for the determination of different characteristics of 
the manifold's dynarnics-without additional physical as- 
sumptions (c.f. RefsOlij). The dypamical FRG starts 
from a Martin-Siggia-Rose actionc3 obtained from the 
equation of motion Eq. (|J) and proceeds with the elim- 
ination of large momentum fluctuations. Thereby the 
parameters entering Eq. (||) are renormalized and char- 
acterize a system for which disorder is less and less rele- 
vant. Finally, the flow is cut off at a scale where the effect 
of disorder can be taken into account perturbatively. 

Our starting point is the system of equations derived in 
Ref. [l5] describing the renormalization of the parameters 
entering Eq. (|4|), to lowest nontrivial order 



$Ai(u) = (e~2C)A l (u)+<:uA' l (u)+T l A / l '(u)+ e^" 8 

Js>0,s'>0 

x [A{'(u)[A,(A,(«' - a)) -Mu + A, (a' - a))] 
-A{(u- V)A{(u + A,«) (5) 
+AJ(A ; ( S ' + *))[£{(« - Xia') - A[(u + X lS )}] , 

d l ln\ l = 2-(- f eSAffae), (6) 

Js>0 



(7) 

(8) 



d l fi=(2-C)fi+cA 2 e" s A;(A /S ), 

Js>0 

d l lnf l = -9+ I e- s \isA'{'{\is), 

Js>Q 

with Ai(ti) = (A d A d -y c 2 ) Ai(u), X, = (vv)i/cA 2 , ft = 
AdA d ~ 2 Ti/c, and // = // — (r]v)i. The exponents ( and 
8 = d—2+2( describe the scaling of the roughness and the 
energy, respectively. A d is the surface of the unit sphere 
in d dimensions and A denotes the short-scale cut-off of 
the theory. Note the important effect of the dynamics in 
rendering the equations non-local on the scale A; propor- 
tional to the center of mass velocity v of the manifold. 
The main goal of this section is to investigate these equa- 
tions in the limit when the external force acting on the 
manifold is slightly below the threshold force, / < f c . 

We first analyze the system of equations (||) - (||) for 
the case of an infinitesimal velocity v = 0+ and concen- 
trate on low temperatures 
to the static FRG equations 



}S. (||) and (||) then reduce 



DtMu) = (e - 2C)Mu) + CuA[{u) + t}a;» 

+A> l '(u)(A l (0)~A l (u))-Al(u) 2 , (9) 

(10) 



-6T,. 



The flow takes the correlator A; through a special 
point l c a (l/e)ln[e/(3|A£(0)|)], L c a A"V° w 
{ec 2 £ 2 /[j4 ( 2Ao(0)]} 1 /( 4-<1 ), where it becomes singular at 
the origin in the limit T —> 0; this is easily seen from the 
equation 



BiA\'(p) a eA{'(Q) - 3A;'(0) 2 



(11) 



satisfied by the second derivative of the correlator at low 
temperatures; exploiting the fact that Aq(0) < we have 
|Af(0)| « |A{f(0)|e d /[l-(3|A{f(0)|/e)(e d -l)]. The cur- 
vature A"(0) diverging at l c marks the occurrence of a 
non-analyticity at the origin which is reflected in the ap- 
pearance of a cusp in Ai > / c at u — 0. Although the initial 
correlator Ao (u) is usually an analytic and even function 
of the coordinate u with vanishing odd derivatives at the 
origin, the function A;>j c has a cusp with a nonzero slope 
Aj(0+) < at the origin when T = 0. Asymptotically, 
Ai(u) approaches a zero temperature 'cuspy' fixed point 
A*(u) describing the disordered phase with a nontrivial 
roughness exponent £. Assuming that the rough shape of 
the fixed point function A*(u) is assumed at the Larkin 
scale l c we can easily find its characteristics: The width 



r «£exp(-CZ c 



(12) 



of A*(u) follows from integrating the second term in (||). 
Comparing terms in (^) at the origin u = we find 
A*' 2 (0+) ~ eA*(0) and combining this with the rela- 
tion |A*'(0+)|r - A*(0) we find the estimates A*(0) ~ 
e^ 2 exp(— 2(l c ) and 
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|A*'(0+)| w e£e 



(13) 



Let us then analyze the force flow (Q) in the light of 
these results. The scale l c divides the flow into two dis- 
tinct regimes, the Larkin regime at small scales I < l c and 
the random manifold regime (/ > l c ). For I < l c we have 
AJ(0) = AJ(0+) = and the force /; obeys the equation 
dji = (2 - Qfi, i.e., fi = e ( 2 -0^/ grows exponentially. 
At the point / = l c the integral term on the right-hand 
side of Eq. (^) jumps from zero to a finite value, since 
the slope of the correlator at the origin does not van- 
ish any longer. If this contribution overcompensates the 
scaling term, i.e., cA 2 |AJ c (0+)| > (2 - C)e (2 ~ C) ' c /, the 
force will start renormalizing to zero while in the oppo- 
site case it will continue to increase. This can be in- 
terpreted in the following way: if the initial value / of 
the external force is smaller than a critical value f c the 
force will eventually renormalize downwards and cannot 
move the manifold. For / > f c the manifold starts mov- 
ing and one should take into account Eq. (^) since the 
problem is not static any more. We therefore come to 
the conclusion that for a system with a dynamics de- 
scribed by Eq. (|J) there exists a finite threshold force f c 
at T = 0. Using the above condition one easily finds 
f c « [cA 2 |A*'(0+)|/(2 - C)]e-( 2 "^ « [£/(2 - C)R/L 2 
with L c = A 1 e lc . Note that the expression for f c 
coincides -with the result obtained from simple scaling 
estimates!! 

At finite temperatures T > it is not possible any 
more to define the critical force density f c as the thres- 
hold below which there is no center of mass motion of 
the manifold, as thermally activated jumps lead to an 
average velocity v > at any finite force /. Of course, at 
low temperatures the velocity is exceedingly small, given 
that it obeys an Arrhcnius-type law. Therefore the time 
needed to observe this velocity might well exceed the time 
scale of the experiment, i.e., from an experimental point 
of view the critical force density still exists with the thres- 
hold f c separating creep-type motion from viscous flow. 
On the mathematical level the nonexistence of the criti- 
cal force density can be explained as follows: at T > 0, 
the slope A'(O-I-) remains zero beyond l c , Aj >; (0+) = 0, 
and the renormalized force density will continue to grow 
beyond the length scale l c even if the initial force density 
/ is smaller than f c . The flow of A; then has to be in- 
cluded into our consideration and the renormalization of 
fi will be found to stop at a larger scale. 

Let us then analyze the flow of the correlator A;, 
Eq. (Q), at finite temperatures in more detail. The non- 
localities introduced by the finite value of A; in (||) can 
be neglected as long as A/ is smaller than the length scale 
introduced by the finite temperature, and we can there- 
fore continue to use the quasi-static equation Eq. (||). 
Below, we will make use of the flow equations only in the 
regime where this condition holds. We also neglect the 
disorder contribution to the temperature renormalization 
in (Q) since it does not influence the main result to the 
accuracy desired here. 



At finite but small T > the correlator flow below 
l c does not differ much from the zero temperature case. 
However, at l c no cusp occurs at the origin — rather, the 
correlator remains rounded on a characteristic scale uj . 
Assuming that outside the thermally dominated region 
close to the origin the correlator has approached its zero 
temperature fixed point shape A* (it) we may estimate 
uj from Eq. (||) by equating the third and fourth terms 
on the RHS, 



T, 



|A*'(0+)| 



A d A d - 2 Te- 



£♦1^-90-1.), (14) 
U r 



with U c the typical elastic energy on the Larkin scale L c , 
U c « (e/Aj) c^ 2 L^~ 2 . Obviously, for low temperatures 
T <C U c , the thermal rounding of the cusp involves a 
scale uj much smaller than the width £* of the correla- 
tor. 
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FIG. 2. Numerical integration of the curvature |A"(0)| 
in the force correlator for d — 1, see Eq. (^|). Four ini- 
tial temperatures To < Td p below thermal depinning at 
Td P = Uc = c£ 2 /L c with values To = T dp /a, a = 40, 20, 10, 5 
have been chosen, while the high temperature curve starts 
with To = 20Td p . The initial growth through the Larkin 
regime involves the exponents e = 3 and e — 5fth = 1/2 at 
low and high temperatures, respectively. Beyond the Larkin 
regime the flow of A"(0) oc — 1/Zj is characterized by the tem- 
perature exponent 9 = 2f — 1 = 0.2495 (c.f., ( |l5| ) with a wap-i 
dering exponent £ given by the one-loop fixed-point valueEd 
(" « 0.2083e). The crossover at the Larkin scale is sharp and 
independent of temperature below Td P ; the uniform vertical 
spacing of the asymptotic curves reflects the temperature in- 
dependence of the energy scale U c . For high temperatures 
above Td p the crossover is shifted to l c (T) 3> l c and the pin- 
ning energy U C (T) depends on temperature, U c ~ T (note 



that I A 



(0)1 



The finite temperature curvature of the correlator fol- 
lows from a comparison of the term TJA" with either the 
fourth or the last term in (||) 
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Af >io (0)«- 



A*'(0+) 2 
ft 



U, 



±J{l-lc). 



T 



(15) 



the finite temperature fluctuations thus regularize the 
'cuspy' divergence occurring at T = 0. The above es- 
timates are correct up to numerical prefactors only; a 
more rigorous derivation can be found in Ref. 

For the following it is crucial to establish that the be- 
havior described by Eq. (051) is valid already soon after 
l c , as the flow equation (g) for A; is quite sensitive to 
the curvature of the correlator at u = 0. Indeed, as 
we have checked numerically (see Fig. ||), after a rapid 
growth within the Larkin regime, the curvature A"(0) 
saturates at a value ~ U c /T with a slow further growth 



due to temperature rescaling. 



Af(0) 



oc T, 



.MI 



As 



long as we are interested in the threshold behavior of 
the barrier close to f c it is sufficient to establish a rapid 
crossover of the curvature from a steep growth below l c 
to a gentle increase above l c (which we will neglect for 
small I — l c > 0). Also note that the crossover occurs 
essentially at the same value of L c independent of the 
temperature T, while the uniform vertical spacing of the 
asymptotic curves reflects the temperature independence 
of the energy scale U c . 

Let us now analyze the system of equations (||) - (||) 
for T > 0, v > close to criticality / < f c . We as- 
sume a fixed deviation f c — f from the threshold and a 
small temperature T <C U c . The velocity then is (expo- 
nentially) small as well and represents the smallest pa- 
rameter in the problem. In the opposite case (T small 
and fixed while f c — f — > 0+) the velocity v oc X 1 /, 8 has 
been argued to scale as a power of temperatureOu, the 
numerical vaLua-of the exponent 5 still being a matter of 
controversyOcIIGa. 

Again, we consider separately the two regimes below 
and above l c . Throughout the Larkin regime the param- 
eter A; remains small and can be set to zero in (§) and 
(7). Furthermore, the temperature dependent term in 
( 3 ) may be neglected initially — it will become relevant 
when A"(0) becomes of the order of U c /T. The flow of A; 
through the Larkin regime then follows from expanding 
the flow equations for A"(0) and A;, (||) and (||), to sec- 
ond order in the small parameter A; and setting TJ = 0, 

ft(ln|Aj'(0)| - 31nAj) =e- 3(2-0 +0(Af). (16) 



Integrating from to l c and using A" (0) 
obtain (up to numerical factors) 



A, 



1/3 



U c /T we 



(17) 



Note that A/ c grows as a power of U c /T as the tempera- 
ture approaches zero, whereas in the following depinning 
regime A; will be exponentially sensitive to T. 

Going beyond the Larkin regime I > l c the function 
A; (u) quickly approaches its fixed point form except for 
a small thermally smoothed region of size uf around the 



origin with the second derivative given by Eq. (||). As 
long as A; < uf one can keep A; = in the integral on 
the RHS of Eq. (||) and a simple integration from l c to I 
provides the result 



A, ~ A ic exp [(2 - C)(i - lc) + ^ (e 9(l ~ h) - l) 



A; c exp 



(I - lc) 



(18) 



where in the last step we have assumed that T/U c <C 

I - l c < 1. 

Turning next to the force equation (^) we note that at 
finite temperature the disorder contribution adds in only 
at a larger scale l^ > l c where A; becomes of the order of 
the thermal rounding uf , in contrast to the zero temper- 
ature case where disorder jumps in at l c . The condition 



\i d cul^f ld /\A*'(0+) 



(19) 



then determines a relation connecting the crossover scale 
Id with the initial velocity v, 



v 



4/3 



exp 



Uc 

T 



(Id - lc) 



(20) 



Below Id the flow of fi is determined by the scaling term 
alone, difi = (2 — ()fi, and a simple integration gives 
/; = /exp[(2 — 0(Ij where we have dropped the small 
correction nv in the definition of /,/ = / — f]v. At Id 
the disorder correction turns on rapidly and we enter the 
depinning regimeLLj at scales I > Id- In this regime we can 
substitute the fixed point correlator A* for A; since now 
A; 3> uj . Furthermore, since still A; <C £* we can set 
A; equal to 0+ in Eq. and we obtain a disorder cor- 
rection cA 2 A*'(0+). As argued before when determining 
the threshold force f c , we have to prevent the force // 
from running away to ±oo and thus the disorder term 
has to match the scaling term; we then arrive at a sec- 
ond relation expressing Id in terms of the applied driving 
force /, 



fe 



{2 _ C)ld ^ cA 2 |A*'(0+)| 
2-C 



= fc^ 2 -^ 



(21) 



With Id close to l c we can expand, ld — l c — (2 — C) 1 (1 — 
/ / /c) *C 1 and combining with ( pQ ) we arrive at the final 
result for the average velocity v 



v oc exp 



Uc fc-f 

T fc 



(22) 



where we have dropped an unessential numerical factor in 
a redefinition of U c . Also, our analysis is not sufficiently 
precise to specify the pref actor. 

Summarizing, we find that close to threshold with 
T/U c -C 1 - f/fc < 2/(1 + fi) the velocity obeys an 
Arrhenius-type law with an energy barrier decreasing lin- 
early on approaching f c . On the other hand, the usual 
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glassy behavior U(f) ~ U c (f c / f)^ is valid at small forces 
f/fc < 2~ 1 ^. In typical experiments the measured bar- 
riers are related to the waiting time t in the experiment, 
U(f) ~ T\n(t/to), and only a limited regime of forces 
with barriers 5 < U(f)/T < 30 is available. This regime 
is, by making use of an extended temperature interval, 
still sufficient to probe both the linear and the glassy 
regimes close to threshold and at low drives, respectively, 
see Fig. § 



glassy 





highT 












low T 











linear 







FIG. 3. Effective creep barrier U at low temperatures as a 
function of external force /. The thin line follows the interpo- 
lation formula U(f) ~ U c [(f c / /) M — 1], properly interpolating 
between the glassy and linear regimes at small drives and close 
to threshold, respectively. The slow relaxation governed by 
the logarithmic decay law U(f) ~ — Tln(l + t/to) limits the 
experimental window to the interval 5 < U/T < 30 depend- 
ing on temperature. Typically, a low (high) temperature mea- 
surement then probes the linear (glassy or non-linear) regime. 



III. CREEP AT HIGH TEMPERATURES 



In this section we consider thermal creep of d = 1, 2 
dimensional elastic interfaces (strings and interfaces mov- 
ing in one transverse direction) at high temperatures 
and for small driving forces. We show that the velocity- 
force characteristic exhibits a non-Arrhenius type behav- 
ior v oc exp [— (f c (T)/f)^]. This dependence derives from 
a creep-type motion with a renormalizcd activation bar- 
rier of the order of temperature, U C (T) ~ T, and involves 
only the renormalized critical force f c (T). For d = 3 the 
parameters f c (T) and U C (T) 3> T depend on tempera- 
ture as well, leading to a non-Arrhenius type creep which 
is non-universal, however, with a resu lt dep ending explic- 
itly on the chosen cutoff. In Section III A we show how 



to calculate the renormalized energy barrier U C (T) and 
the threshold force density f c (T) using scaling arguments 



and then present a mo re rig orou s ana lysis using (dynam- 
ical) FRG in Sections [III E| and [HIC| . 



A. Scaling Analysis 

While the mean squared displacement (it 2 )th ~ TA/irc 
is bounded in d = 3 (with A -1 the intrinsic cutoff of the 
manifold), the thermal displacement (or wandering) of 
strings and interfaces (d = 1, 2) grows unboundedly, ei- 
ther with distance z or time t, 



(u 2 (z,t)) th = ([u(M) -u(0,0)] 2 ) th 



(23) 



2T f d d q 



(2nf 



q~ 2 Re 



I — e *q- z e -( c /'))'? * 



-[z 2 + (c/^] 1/2 , d = l, 

7TC 

In [A 2 (z 2 + {c/r,)t)], d=2- 

the pinning length L C (T) set by the disorder landscape 
then has to provide the necessary cutoff which in d = 3 is 
given by the intrinsic cutoff A -1 (note that the g-integral 
in ( pH ) is dominated by small (large) q for low (high) di- 
mensions; hence, depending on the dimensionality of the 
system the amplitude of thermal fluctuations is deter- 
mined by short (d — 3) or long (d = 1, 2) scales). This 
implies the existence or absence of a separation of scales 
for thermal and disorder effects: While for d = 3 these 
scales are separated, A -1 <C L C (T), no such separation is 
effective in dimensions d — 1, 2; thermal effects smearing 
the disorder landscape are active on scales L < L C (T) 
while disorder takes over for L > L C (T), i.e., at the same 
scale. The temperature induced smoothing of the disor- 
der potential then follows different rules in low and high 
dimensions, as we are going to discuss now. 

In order to discuss pinning and creep we have to de- 
termine the renormalized disorder landscape. Assuming 
pinning to involve longer time scales than thermal fluc- 
tuations we average the pinning potential over thermal 
excursionsEl, 



«£pin(£))t) 



dt 



d a z 



t Ud ±\d* z 

to Jo to J 

x{V(u(z,t),z)V{u(z',t'),z')} 

U 'r f° ^K[Ht)-u{t')\] 

o r o Jo T o 

( f 2 \ 1/2 

LdK (°) ,..2,K-1 > ( 24 ) 



L 



(m 2 (A _1 , L))th 



where the mean squared thermal fluctuations 
(tt 2 (A _1 ,i))th are cutoff by A" 1 or L = L C (T) in high 
and low dimensions, respectively. Here, K(u) denotes 
the potential correlator which is related to the force 
correlator A(u) used above via —K"(u) = A (it). The 
result (|3) tells us that at high temperatures, thermal 



G 



fluctuations replace the basic length scale £ of the dis- 
order landscape by the scale (it 2 (A -1 , L))\^ > £ (note 
that the energy scale of the disorder potential remains 
unchanged). Comparing this smoothed pinning energy 
with the elastic energy c(u 2 (A _1 , L)) t hL d ~ 2 ~Twe ob- 
tain the new pinning scale replacing the T = Larkin 
length L c = (c 2 £ 4 /^(0)) 1/(4 ~ d) , 



L C (T) ~ L c 



V(u 2 (A-\m 5/2 ' 



l/(4-d) 



Km 



(25) 



For high dimensions we find L C (T) ~ L C (T/T dp ) 1 / 2 ^ F 
with the Flory exponent £ F = (4 — rf)/5 and the depin- 
ning temperature Td p = c£ 2 /A (this result follows from 
simple scaling u oc and u oc T 1 / 2 ). In dimensions 
d = 1, (2) the corresponding result takes the form (here 
we concentrate on the case n = 1 , see Ref . ^9] for a dis- 
cussion of the marginal situation in d + n = 1 + 2) 



L C {T) ~ L c 



T 



(26) 



with the temperature exponent A = 5, (5/4) and the de- 
pinning temperature Td p = (cif (O)^ 2 ) 1 / 3 , (c£ 2 ). Com- 
paring with the above Flory exponent wc sec that thermal 
fluctuations indeed are much more important in d = 1, 
while for d — 2 the corrections are only logarithmic (not 
shown in (p6|)). 

The energy barrier and the threshold force are renor- 
malized correspondingly; for d — 3 we have U C (T) ~ 
U C {T/T dp ) 7 / 2 » T, with [7 C = d 2 L c and T dp = c£ 2 /A, 
while for d = 1, 2 the barrier 'saturates' above Td p , 
U C (T) ~ T. The critical force density is renormalized 
according to f c (T) ~ / c (T dp /T) 9 / 2 in d = 3 and takes 
the form 



/ccn 



(^(^c^)))^ 2 

L C (T) 2 



fc ( y 



(27) 



with the temperature exponent k = 7, (2) in dimensions 
d = 1, (2), restricting ourselves again to the case n = 1. 
Using the usual creep formula U(f)/T = {U c /T){f c J fY 
and inserting the new temperature dependent values for 
U c — ► T and / c — ► f c {T) we obtain the creep exponent 
(fc/ f) fJ '{Tdp/T) KfJ ' producing a non-Arrhenius type creep 
in d = 1, 2 involving only the renormalized critical force 
density f c {T). For higher dimensions d > 2 the renormal- 
ized barrier U C (T) 3> T remains large, a consequence of 
the separation of scales mentioned above. In the follow- 
ing we rederive these scaling results via the more rigorous 
analysis provided by the (dynamical) FRG scheme. 



B. Functional Renormalization Group 

In a first step we rederive the crossover scale L C (T) 
via the functional renormalization group. The nonlinear 



terms in the flow equations are still small during the ini- 
tial stage of the RG flow — neglecting them we first solve 
the linear equation. The length L C (T) then appears as 
the characteristic length where the nonlinear corrections 
become of the order of the linear terms. The analysis 
is conveniently carried out for the potential correlator 
Ki(u) which follows the flow equation 



+ \K{'{u) 2 



diKi(u) = (e - K)Ki{u) + QuK[(u) + f,£,"(u) 

k['{u)k['(Q), (28) 

while the temperature flow is given by 

am = -e%. (29) 

Using the ansatz 

Ki(u) = exp [(e - 4C - (2 - d)/2)l] hiue 81 ' 2 ) (30) 

the linear part of the flow transforms into a Fokker- 
Planck equation describing the probability distribution 
Pi(u) for an overdamped particle moving in a parabolic 
potential at constant temperature TqEj, 



diPi(u) 



-d u [uPi{u)]+f dlPi(u), 



(31) 



for which the fundamental solution is well-known. Solv- 
ing the initial value problem for Pi and inserting in Eq. 
(pp|), we obtain 



Ki(u) 



exp [(4 - d - 5QI] 
du' Kq(u') exp 



2-d 



4wT e- el (l-e-( 2 - d ) 1 ) 
(2-d)(u-u'e-«) 2 



1/2 



4T e- m (l-e-( 2 - d ) 1 ) 



For high temperatures such that f Q e~ el (l-e~ ( - 2 ~ d)l )/(2- 
d) > ^ 2 e~ 2 ^ 1 the factor ko(u') acts as a 5-function and 
can be extracted from the integral, 



k(u) 



= (3-5C)i 



exp 



4T e- ei 



J Kq (u) du 



,(2-50* 



(inToe- 61 ) 1 / 2 
exp 



4T le- ei 



,(1-5C)* 



(47r7ye- e 1/2 
exp 



d= 1, 



d= 2, 



4T e- 2 « 



(47rT e- 2< ; / ) 1/2 



(32) 



The above constraint on the temperature simplifies to 
f > e"^ 2 , f > £ 2 /l, and f > £ 2 for d = 1, 2, 3, re- 
spectively, where £ is the initial width of the correlator 
Kq; expressing To and I through the physical quantities 
T and L we recover the condition (u 2 (A _1 , L)) t h > S, 2 ■ 
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The solutions (|32|) still involve the wandering exponent 
£. Although the final physical results do not depend on 
the particular choice, it is a matter of convenience to 
adopt the thermal values £th = 1/2, t h = for d = 1, 
Cth = 0, th = for d = 2, and Cth = 0, 8 th = 1 for 
d = 3 and have the correlator flow towards a thermal 
fixed point. The renormalized correlator ([32]) then be- 
haves very differently for large and small dimensions: In 
d = 3 the transverse scale u does not change (as we chose 
C = 0) and the initial correlator of width £ is replaced 
with a new correlator of width (u 2 ) t h ~ To = TA/c. 
This contrasts with the situation in d = 1 where u does 
rescale (as we chose £ = 1/2) and the physical width of 
the correlator increases with I to follow the mean thermal 
displacement amplitude (u 2 (L))th ~ TL/c. For d = 2 the 
physical width grows only logarithmically. 

The flow (||^) indicates that the thermal fixed point 
is unstable as the amplitude of the disorder grows ex- 
ponentially under the FRG transformation. As the 
nonlinear terms in ( |28| ) become large beyond the scale 
L C (T) — A _1 e' c ( T ^ we cannot neglect them any longer 
and the flow crosses over to approach the disorder dom- 
inated fixed point (the wandering exponent £ then has 
to be modified accordingly, ( = 0.2083e for random bond 
disordei£3). The pinning length L C (T) replaces the T = 
Larkin length L c and can be found from a comparison of 



linear and quadratic terms in the flow equation (28), 



iv ;c(T) (0)^ (T) (0) 2 



(33) 



making use of the result (p2|) . It is easily verified that 
the crossover condition (|35|) together with the explicit 
solution ( B2[) o f the linearized flow equations then yields 
the results(|25|) for the crossover length L C (T) obtained 
above with the help of scaling arguments. Note that for 
d = 1, 2 the linear term TiK" in (|2|) gains in impor- 
tance as we integrate through the Larkin regime, hence 
thermal rounding persists on all scales I < l c (T). On the 
contrary, for d = 3 the thermal rounding term is most 
important at small scales I ~ 1 where it quickly replaces 
the width £ of the correlator by the mean thermal dis- 
placement amplitude (u 2 )^ 2 ; upon further scaling the 
effective temperature decreases and at the Larkin scale 
the thermal term is down by a factor of e~ la< ^ T K This 
again reflects the different role the temperature plays for 
different internal dimensions of the elastic manifold. 



C. Dynamic Functional Renormalization Group 

Finally, let us see how the high temperature creep (as 
characterized by the temperature dependent critical force 
density f c (T) for d= 1, 2) appears directly from the dy- 
namic FRG treatment. The analysis parallels the treat- 
ment at low temperatures, however, we have to be more 
careful in distinguishing between the cases d = 1, 2 and 
d = 3. 



Let us analyze the flow of the force correlator A/(u) = 
-k['{u). Following the full flow up to l c (T) the force 
correlator assumes a shape with a height and width as 
given by (|32|) . At l c (T) the non- linear terms in the flow 
equation ( p8| ) have become important; beyond l c (T) the 
correlator quickly flows towards the disorder dominated 
fixed point function A* (it), the linear temperature term 
Ti A/(u) smoothing the flow in a region of size uf around 
the origin. Assuming again that the fixed point func- 
tion A* (it) derives is rough shape from the correlator 
Ai c (T)( u ) a t crossover, we can use the result ( [32] ) in com- 
bination with the flow equation (Q) to find the character- 
istic features A*(0), £*, and A*'(0+) of the fixed point 
function and the rounding parameters AJ'(O) and uf of 
the cusp. 

Using (|3^) and the crossover condition ( |33| ) we find 
the height A*(0) = -i^ (T) (0) « f . The slope of 
the fixed point function at u = 0+ again follows from 
comparing terms in the flow equation, |A*'(0+)| w 
A*(0y/ 2 « (f ) 1/2 and we find the width £* w (To) 1 / 2 

w (u 2 (L c {T)))H 2 e-^ MT \ see also the result @. 
The width uf of thermal rounding derives from U; ~ 
f//|A*'(0+)| and we find the result uf /£* - e -e(i-UT)) 
in d = 1, 2, with 8 = d — 2 + 2£ and £ the random mani- 
fold exponent; thus at l c (T) the width of thermal round- 
ing equals the width of the correlator, uf w £*. The 
curvature A"(0) ~ A*(0)/T) is correspondingly small, 

- P (l-lc(T)) 



A"(0) ~ —eU-'etJ;^ see Yig. 2; comparing this result 
with Eq. (15) we conclude that the barriers 'saturate' to 
follow the temperature, U C (T) w T. This is quite dif- 
ferent from d — 3: Here, the thermal rounding affects 



only the narrow regime u 



-l c (T)„-9(l-l c (T)) 



(U C (T)/T) e eiyl Z <=( T )) around the origin and the curvature 
is already large at l c (T), A/(0) » - e ^ T ) e e ('-^ T » = 
-(£/ c (T)/T)e e «-^ T », where U C (T) = c(u 2 ) th L c (T) » 
T. Thus, as anticipated above, the curvature of the cor- 
relator at l c (T) is thermally reduced to the order of 1 in 
d = 1, 2 while it is large in d = 3. Hence for d = 3 the 
situation at high temperatures T > Tdp is not different 
from that at low temperatures. 

Next we integrate the flow for the velocity and force 
parameters A/ and //. We determine the creep scale ld(T) 
twice, using the condition A; d (y) = u f d (T) ^° re l a t e the ve- 
locity v and the scale ld(T) and a second time from the 
onset of the disorder term in the force equation ( pl| ) , pro- 
viding a relation between / and ld(T); combining these 
results we obtain the desired velocity-force characteris- 
tics. In doing so, we have to be careful to use the above 
high temperature estimates for A*'(0+) in Eqs. (|l9| ) and 

Integrating the flow equation @ for A; through the 
Larkin regime and then up to ld{T) we find the first re- 
lation 



v cx exp[-e e (W T )-^ T »] 
Integrating next the force equation m 



(34) 



we obtain 
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fi d (T) ~ fe 



(2-( th )l c (T) + (2-C)[l d (T)-l c (T)] 



(35) 



and equating this to the disorder induced term 
cA 2 |A*'(0+)| in the flow equation (||) we arrive at the 
second relation 



/ ~ / c (T)e- (2 ~ c)MT) ^ (T)1 



(36) 



with the critical force density / C (T) = cT^ 2 A 2 
e -(a-Cth)I.(T) „ c(u 2 (L c (T)))H 2 /L 2 c (T}, in agreement 



with fl27|). Combining the results ( p4| ) and (|36| ) we 
find the velocity-force characteristic describing the non- 
Arrhenius type creep at high temperature, 



exp 



V oc < 



exp 



fjry\ e/i2 -° 

0/(2-0 



1, 2, 



d = 3: 



In conclusion, using dynamical functional renormaliza- 
tion group theory we have derived the linear scaling of the 
creep barriers close to f c and have put the non-Arrhenius 
type high temperature creep of low-dimensional mani- 
folds on a firm basis. The simple behavior of the creep 
barrier close to threshold appears surprising — consid- 
ering the 'non-trivial' threshold exponents due to a di- 
verging nucleus obtained for elastic manifolds trapped in 
a washboard potential (see Ref. [|) one is tempted to ex- 
pect a non-trivial exponent for the random case as well. 
However, from our analysis we conclude that there is no 
new diverging scale associated with creep near threshold. 
The linear decay of the creep barrier then follows from 
a regular expansion and no effects of critical fluctuations 
are picked up. 

We thank Pascal Chauve for discussions and the Swiss 
National Foundation for financial support. 
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